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For almost twenty years the effects of a nonzero strange quark mass on the equation 
of state of cold and dense QCD were considered to be negligible, thereby yielding only 
| minor corrections to the mass-radius diagram of compact stars. By computing the ther- 
modynamic potential to first order in a s , and including the effects of the renormalization 
group running of the coupling and strange quark mass, we show that corrections can be 
of the order of 25%, and dramatically affect the structure of compact stars. 

in ■ 
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Compact stars provide the most promising "laboratory" to constrain the equation of 
state for strong interactions. In this work we compute the thermodynamic potential for 
cold quark matter with two light (massless) flavors, corresponding to the up and down 
quarks, and one massive flavor, corresponding to the strange quark, in perturbation theory 
to first order in a s in the Ms scheme [[Q. In this fashion, we can easily include modern 
renormalization group running effects for both the coupling and the strange quark mass. 
^ ■ We find that the corrections due to the running nonzero mass are sizable, and should 
not be neglected in the evaluation of thermodynamic quantities. Solving the Tolman- 
Oppenheimer-Volkov (TOV) equations using our equation of state in the presence of 
electrons, we show that a running strange quark mass dramatically modifies the mass- 
radius diagram for quark stars, even at first order in ot s . 

The thermodynamic potential for cold QCD in perturbation theory to ~ a 2 was first 
computed a long time ago [El Bill!- Nevertheless, the original approach to quark stars [ 
03 IE] made use of the bag model with corrections ~ a s from perturbative QCD to compute 
the thermodynamic potential. In the massless case, first-order corrections cancel out in 
the equation of state, so that one ends up with a free gas of quarks modified only by a 
bag constant. Finite quark mass effects were then estimated to modify the equation of 
state by less than 5% and were essentially ignored. 

A few years ago, corrections ~ a 2 with a modern definition of the running coupling 
constant were used to model the non-ideality in the equation of state for cold, dense QCD 
with massless quarks [IHUHl- This approach can be compared to treatments that resort to 
resummation methods and quasiparticle model descriptions [EH HUES EH]- Remarkably, 
these different frameworks seem to agree reasonably well for fi » 1 GeV, and point in 
the same direction even for fi ~ 1 GeV and smaller, pushing perturbative QCD towards 
its border of applicability. 

Even the most recent QCD approaches mentioned above generally neglected quark 



2 



Eduardo S. Fraga 



masses and the presence of a color superconducting gap as compared to the typical scale for 
the chemical potential in the interior of compact stars, ~ 400 MeV and higher. However, 
it was recently argued that both effects should matter in the lower-density sector of the 
equation of state [ I14j. In fact, although quarks are essentially massless in the core of 
quark stars, the mass of the strange quark runs up, and becomes comparable to the typical 
scale for the chemical potential, as one approaches the surface of the star. 

In what follows, we present an exploratory analysis of the effects of a finite mass for the 
strange quark on the equation of state for perturbative QCD at high density, leaving the 
inclusion of a color superconducting gap in this framework for future investigations. To 
illustrate the effects and study the modifications in the structure of quark stars, we focus 
on the simpler case of first-order corrections. Results including full corrections ~ a 2 , 
as well as technical details of the calculation at each order and renormalization, will be 
presented in a longer publication [IT5]. 

The leading-order piece of the thermodynamic potential of QCD for one massive flavor 
is given by [0E1H 
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where N c is the number of colors and u = y/ fi 2 — m 2 . 

Using standard quantum field theoretical methods, one obtains the complete renormal- 
ized exchange energy for a massive quark in the Ms scheme (in the limit T — > 0): 
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The thermodynamic potential to order a s for one massive flavor, given by the sum of 
Eqs. (J21) and (0), depends on the quark chemical potential fi and on the renormalization 
subtraction point A both explicitly and implicitly through the scale dependence of the 
strong coupling constant a s (A) and the mass m(A). The scale dependencies of both a s 
and m, which in the following we will take to be the mass of the strange quark, are known 
up to 4- loop order in the Ms scheme [El- Since we have only determined the free energy 
to first order in a«, we choose 
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where L = 21n(A/A M g), (3 = 11 
The scale A 



2JV//3, and ft = 51 - 19iV//3 and we take N f = 3. 
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Mg and the invariant mass m s are fixed by requiring [ El ol 



100 MeV at A = 2 GeV; one obtains A- 
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380 MeV and m s ~ 262 MeV. With 
these conventions, the only freedom left is the choice of A. 

To study the effect of the finite strange quark mass on the equation of state for electri- 
cally neutral quark matter with 2 light (massless) flavors (up and down quarks) and one 
massive flavor (strange quark), we have to include electrons, with chemical potential \x e 
and assume beta equilibrium. Since neutrinos are lost rather quickly, one may set their 
chemical potential to zero, so that chemical equilibrium yields fid = Us = A* an d + A^e = 
[i , with fi u , fid and fj, s the up, down and strange quark chemical potentials, respectively. 
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On the other hand, overall charge neutrality requires (2/3)n n — (l/3)n<2— (l/3)n s — n e = , 
where is the number density of the particle species i. Together, the above equations 
insure that there is only one independent chemical potential, which we take to be fi. 
Number densities are determined from the thermodynamic potential by rii = —(dQ/dfi) 
and the total energy density is given by e = Q + J2i Hi^i, where = J2i (jn + ^i^) 
and again i refers to the particle species. The pressure is P = ns-£^ — e , where 
n B — | { n u + n d + n a ) is the baryon number density and the Gibbs potential per par- 
ticle is given by = (u. u + u, d + u. s ). We restrict the freedom of choice for A(/z u , u, d , u, s ) 
by requiring that in case of vanishing strange quark mass all quark chemical potentials 
and number densities become equal so that p( ms= °) = — ^( m s=°) and, consequently, one 
has /i e — > 0. Furthermore, in order to compare our findings to existing results in the 
literature [El El, we require that in the massless case A = 2fi. Consequently, we choose 
A = | + fid + fi s ), but have tested that our results are not much affected by other 
choices obeying the above conditions. 

The effects of the finite strange quark mass on the total pressure and energy density 
for electrically neutral quark matter (plus electrons) are given in Fig. ^ There we show 
results for 3 light flavors and running coupling, corresponding to the case considered in [ 
ISj, and for 2 light flavors and one massive flavor, with both running coupling and strange 
quark mass (which reaches m s ~ 137 MeV at fi = 500 MeV). 
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Figure 1. Pressure and energy density scaled by Stefan-Boltzmann values for A = 
| (a^u + fid + l^s)- We show results without renormalization group improvements (dash- 
dotted lines), running coupling (dashed lines), both for m s = 0, and results with running 
mass and coupling (full lines). 

As can be seen from this Figure, there is a sizable difference between zero and finite 
strange quark mass pressure and energy density for the values of the chemical potential 
in the region that is relevant for the physics of compact stars. As has been noticed by 
several authors [HHTUHE], the resulting equation of state, e = e(P), can be approximated 
by a non-ideal bag model form e = 4:B e ff + aP. Here a ~ 3 is a dimensionless coefficient 
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while B e ff is the effective bag constant of the vacuum. Concentrating on the low-density 
part of the equation of state, one finds for massless strange quarks the parameters -B^// — 
117 MeV and a ~ 2.81 while the inclusion of the running mass raises these values to 
-Be// ~ 137 MeV and a ~ 3.17 (all values having been obtained by including a running 
a s in the equations of state). Therefore, we expect important consequences in the mass- 
radius relation of quark stars due to the inclusion of a finite mass for the strange quark. 

The structure of a quark star is determined by the solution of the TOV equations. 
Corrections to the mass and radius of quark stars due to a running strange quark mass 
can be very large, ~ 25% [H]. 

Also, while the most massive star for the m s = equation of state (with M/M & ~ 3.2 
and radius ~ 17 km) has a central density of % — 0.5 fm -3 , this number increases to 
n B ~ 0.83 far 3 (at fj, = 470 MeV) for the heaviest star (M/M ~ 2.16 at ~ 12 km) of the 
massive equation of state. The inclusion of ~ a 2 s corrections to the pressure will increase 
its non-ideality and produce quark stars which are smaller, denser and less massive [ 
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